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ABSTRACT 


The behavior of the system hazard rate function of a two 
Component parallel system is investigated. The inter- 
relationships between the probabilities that the components 
composing the system are alive and the system hazard rate is 
examined with special attention to certain points where 
there are important changes in the behavior of the hazard 
rate function. The behavior of the system hazard rate 
function is shown to depend upon the rates of change of the 


probabilities that the components of the system are alive. 
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i. v DNSRODUCTION 


in thie itele Of rettarcaruy, co nere has been an extensive 
study of the properties of the increasing hazard rate (IHR) 
and decreasing hazard rate (DHR) classes of life distribu- 
mloms., tne expomenvial tire  adisuribuvions with their 
constant hazard rates form a dividing line between these two 
classes. It is reasonable to expect that the life distribu- 
tion of redundant systems composed of devices with exponen- 
Pewee ite Giseribueitons fall into the Lik class. However, 
when two devices with different exponential life distributions 
Pew conoiicd inea waralleole reG@undanemsyistelml .Glc system 
Mezarad rate is at first increasing and then decreasing. The 
underlying causes of this interesting behavior is the subject 
oe tne tChnesis. 

The eee qderineleom, ol ume Drobapltmns tele naZzard skate 
maneulon wot a device with a random litetime LT willbe used 
throughout this thesis. The hazard rate r(t) of a device 
is the density ranee ton f(t) for the lifetime of the device 
divided by es survival funce lene er =F ioe |. 

A Seinen has an increasing hazard rate if r(t) is 
monotone non-decreasing in t, and has a decreasing hazard 
rate if r(t) is monotone non-increasing int. References 1 
and 2 contain an extensive survey of the properties of the 
THR and DHR classes of life distributions. The exponential 


life distributions are included in both of these classes. 


} 





mn Reference 3 ,e8SaryeanGd -eoseman Save Suiticremt conditions 
mom a SYStEM CO Mavewan Jmerecasing hazard rate when it is 
ceompescdmot, 1denbiealacomponentcawivamamcreasing hazard 
rates. A parallel system of two components, with identical 
exponential life distributions, satisfies these conditions 
and it has an increasing hazard rate. Esary and Proschan, 
also gave an example to show that when two components with 
different exponential life distributions form a parallel 
system, the system survival function need not have an 
increasing hazard rate. In this case the system hazard rate 
initially increases, overshoots and then decreases to an 
asymptotic value equal bo the lowest of the two component 
moire rates. 

ticoeccaon Il, the derivawron of the system hazard rate 
for a system of two parallel components having constant 
non-identical failure rates is reviewed. In Section III, 
a characterization of the system hazard rate as a function 
Or some conditional state probabilities is presente and 
their interrelationships are discussed. In Section IV, an 
empirical approximation for the time at which the system 


hazard rate is a maximum is exhibited. 





IT. SYSTEM HAZARD RATE 


Consider a system of two independently functioning 
components with exponential lifetimes Ti sb and survival 


muineoLOns 
F(t) ae 
F(t) a 


Assume without loss of generality that A5?A,- The equations 
for the system survival function and the system hazard rate 
are generally known. 


Math the components in parallel, the system lifetime T 


has a survival function 


F(t) = P[max T, »T5°t] 
7 “A,t Az —(A, +A, )t 
=e +e - =e 
The density of T is 
_ qd 
1 ae Go) | | 
-A.t Sy (ea Me 
- 1 2 ae 
= die tine —(A, +A, )e 
The system hazard rate is 
5 )) = a) 
HiGe) 
~r¢ -r¢ ~(XA. +A, )¢ 
a 2 eZ 
_ Ae tire —(A, +A,)e 
F(t) 





Examining or (t sate t=Usancduas tuaperoOccicsaimiimaty ¢1ves 


mee expected resulgs. 
r(0) = 0 


lim r(t) = min(A, ,A5) = 


£00 


alt 
The derivative of r(t) with respect to time is 


F(t)f'(t) + £°(t) 


1 Gr ee= 
F* (t) 


9 “Ait 5 “Ant —(A 


3 +r, )t 
Lr. e +r e mich, — Age) Je 


ol 
F(t) 


Looking at the sign of r'(t), it is clear that the 


Hazard rate is increasing if and only if 
-rA_¢ -rA.¢ 


Ne 7 i‘e = s (A, -Ay)® eld decreasing i rand Only sit 


-r,¢ -A,t 


fe 7 thie 2 < (A,-A,)® . If AQHA, » v(t) 4s 


2 & a ia 


increasing for all t>0. The system hazard rate must be at 
a maximum for the value of t which is a solution to 
~A.¢ -A¢t 


2 al 2 2 i 
do e +h e ~(A,-A5) = 0 





TET. SMS teVOHAZARD RATE 
UN TERMS. OF 
Swe Phnom As Tha aS 
If the parallel system is alive at some time t70, it 
exists in one of three possible states, i.e. both components 
mye, COMPONCM OMe alivewane cOMmponeny two dead, or 
component one dead and two alive. The system hazard rate 
and its derivative Paver summrececxpresstOns in terms of the 
conditional probabilities for these states and their 
ferivatives. The probability that the system is in state 
ieea=O,1,2¢, given it is alive at a time t, will be called 
the state SuObaD Lee, denoted by P,(t). For a system with 
lifetime T and component lives T) and Ts; the state 


probabilities are: 


(1) both components alive 


Po(t) = P[T,>t,T,>t|Tot] 


P[T, >t ,T,>t,max(T, ,T,>t) | 
——Fimax(l»Tp5b)1 
PLmax Bj sts>t 


—CA, tA, )t 
e€ 


Gey) 
(2) component 1 alive and 2 dead 


(te) SPI SG piace Sb! 


P[T, >t ,T,<t,max(T, ,T,>t) J 
P[max(T, ,T,>t) ] 





-A,¢ —(A, tA, )t 





1 
P,(t) = == 
Pty) 
sgh 
= ] - 
F(t ) 
ee 


iH 


Ce : -1)P)(t) 3 


(3) component 2 alive and 1 dead 


iH 


Bata) SIC ie Sete aes | 


P{T,<t ,T >t smax(T, ,T,>t) J 


2] 
P{max(T, ,T,>t) J 
“Ant -CA,tA5)t 


= a 


F(t ) 
-r.¢ 
aie 
inGe), 





1 + 


ee 
(e = 


-1)P,(t) 


fies Clear that Pott) + P(t) + P(t) eT ifore elimi c 
peaminane the state probabilities at t = 0 and as t 


approaches infinity produces the following results. 


P,(0) = 1 Jim Pace) = 0 
0 3 ae 0 
eG) = fe lovie2> (Ce eal 
al oes al 
P,(0) = 0, lim Solie ues 
Foo 


10 








At first glance, it might seem that P,(t) would equal 
P(t) at the méan tirime to failure for component 2, hut, 


P(t) = P(t) ais enol weyallliy aLi% 


-A,t —(A, +A, )t —CA, tA, )t 
e -e og FS 
F(t) Go 
Gx 
=a —~(A,+A,)t 
memeeer se ; 
-A~¢ 
or ia 2c - : 
or t = == gaye 2 
do 


ioe probability that a particular single component is alive 
at some time t is greater than the probability that both are 
still alive for any time t greater than the mean time to 
failure of the other component multiplied by a constant ln ée. 


Further, for A >A,» P,(t) > P(t) HOre cullen Sates 


2 
P, (t) > P(t) ii ation y edt 


“Ant “Ajt 
1- = Seelee = 
ECG) ce 








ial 





oT 


OR 


The derivatives of the state probabilities with respect 


to time are: 


—(A, +A,)t ee = 


: F(t)(-(A,+A,)e C0 )) 
EQh eS 
F(t) 
—(A, tA, )t 
= $e Ae EL ata) 1 
F(t) F(t) 


Po(t)ir(t) - (A, #4) 








_ “Ant “rot _. 
7” F(t)(-A,e © ) -e F (t) 
PCs) pte oe 
: F°(t) 
Ant 
= =——— [A, - r(t)] 
F(t) 
ot ae 
be F(t)(-Aje 7 )-e ~ ¥ (6) 
Pete 
2 Fee 
Aye 
= CA, al RiGha)y| 
F(t) 
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P(t) decreases in time from Py (0) = 1 to lim Py (t) = 0, and 


Po(t)<0 for all t. P(t) is increasing i from 
P, (0) = 0 to mi P,(t) = 1, and P(t)>0 P@je ea leet. 

P(t) has a maximum at a time t which can be found by 
etestituting for r(t) in P(t) and setting this expression 


equal to zero, i.e. when 





~(\,+A,)t ~(2).+i.)t 
Cee. 2 web ae os 
ie 2 2 a, 
F(t) 
(A, +A) t 
Peo plying both sides by e and combining terms 
gives 
“r,t - AomAy 
: yk 
2 
and 
Pa 
1 Zi 


The value of P,(t) at maximum is 


7 
ee 
1.2 Seen ee) “I 
9 ~ Aro o7y 
ee Ors Oe ccm eral 6 OE ee 
ee LAS Oe) > ee of 4 


2 I 2 2 


and this is less than one for all hoy>hy 


As a consequence of the definition of the state 


probabilities and the requirement that 


Po(t) + P(t) + Po(t) =1 , 


nS 





the derivatives are related and 
1 1 1 
P(t) P, (t) + P(t) = Q 


Therefore , BC) =~ -[P,(t) + P,(t) and P(t) increases 
when -P,(t) > P, (t) and it decreases when -P,(t) < Parca) 
P(t) increases as long as P(t) decreases in time faster 
than P, (t) increases, reaches its maximum value when 
~P,(t) = P,(t), and then decreases because P(t) is increas- 
ing in time faster than P(t) is decreasing. 

A hazard rate may be thought of as a conditional 
instantaneous probability of failure at time t, given 
SUrVival to time t. The hazard rate for the system, given 
Mm@ee it is in one of its @hreem@etates, is 0, his and 5 
respectively. The system hazard rate may be written in 


terms of these hazard rates and the state probabilities as 


r(t) 


0-Py(t) + AP, (t) 7 A5P5(t) 


ale 2 
=O- + = OF 
F(t) F(t) 


—(A, +A,)t —A,t —(A, +A, )t —r gs ir 


F(t) 


which reduces to the expression previously given for r(t). 
The derivative of the system hazard rate with respect to 


time may be written as 


D Ce) = dae, We) Ge) 
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men the condition for the Maximum point of Chew nazard rate. 


r(t) = 0, requires 
t t 
A P(t) + AoP5(t) = 0 
t t 
S@bstituting for Pi (t) and P5(t) yields 


-Aot -d,¢ 
[An - r(t)} + ay 
iG) 


& 


r 
Ei cee 








D Ee eiGcua), = 30 


Scomvime this for r(t) gives its maximum value 





—-rA.t —-rA.t 
r.r.[e + +e * ] 
Pease 
CLC) ga ae) en =e 
2 1 
Ae + Are 
“Ant | 
$ 
From P.(t) = © Dy SmGeN) = © for aid, che aig elleas: 
1 F(t) 


iets rit) < do Ora et i it selena rane mere (ue) nee 


fond Only if 


-rA,¢ -rA~t 
Noh <> ee 0G 5 
12 >) 
“A56 “A,t il e 
A, e + A5e 
or 
-A.¢ -A~t -Av\~t -r.¢ 
il 2 Zs 2 iL 
Ay Ave tA, A5e > hy é thi Aje ; 


db) 





OR te 


or Ae, 


fimerefore, at its maximum, A, < r(t) < A 


a e 


~A,t 





t e 
From P(t) = 
F(t 


LA, = CC) iiines lt 1S%cleart that 


r 
ao thestimens rt — 2a ( e ) when P(t) reaches 
ak Ay~ Ay 2 


mc) = dX 
Ay 





mecemaximum value. If t > t , then r(t) > dy: Ob.) a> 4, if 


Emad only if 


dP, (t) iF AoPQ(t) > hy 


or AnPo(t) > ae: - P,(t)) : 
or (A5-A, )P5(t) > A P(t) 


since 
| A,t 
P(t) = (e = cc) ae 


tmen r(t) > ry diana ce ‘Onilayar: f 


dit 
(A,-A, ) Ce ~1) >a, , 
a 
amd for C92 1 = _ in — ) 
iL ew i) 


iLie 





Lipemollows “that 


i ae )] 
exes eH TL 
(Ag-Ay, Ce 2S ees (2g=A, ibe 1dAy eo) 4) 
5 
(A5-A,. ) (-———  - 1) 
aa Oa 15 
= hy 


As 5 approaches aT in value t becomes very large, and if 


h, = Ano =X, the system hazard rate r(t) approaches i 


a 
asymptotically from below. For 5 ee hy > tT becomes very 
small and the system hazard rate is greater than hy almost 
miescantly. 

From r'(t) = AyP, (t) + eae) a. PCE) Wee anereascane ut 
and only if AaP,(t) + oP) > 0 , and decreasing if and 
only if AP, (t) + oP) < 0. The behavior of the system 
hazard rate is directly dependent upon the rates of change 
of the state probabilities P, (t) and Po(t). r(t) increases 
as long as P,(t) and P(t) increase and it continues to 
increase as long as AyP, (t) > [oka . For 
me, (t) < [A5P (4) , r(t) decreases. 

Figure 1 shows the state probabilities and the system 
hazard rate plotted as a function of time. As the failure 
rate of component two increases in relation to the failure 
rate of component one, these curves are "pushed towards the 


left", the system hazard rate increases rapidly, and its 


17 





maximum value increases. P(t) Geereasese VOmecCroe Very 
rapidly while P,(t) increases to one almost instantly. 


P(t) becomes more peaked and it reaches its maximum value 


SoOner. 
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IV. EMPIRICAL RELATIONSHIP 
BETWEEN TWO MAXIMA 


An empirical investigation of the relationships between 
the state probabilities and the system hazard rate was also 
undertaken. The system equations were solved for various 
Values of their parameters which permitted comparison of 
Mees resulting curves and verification of the analytic 
memos Of Section II. 


The time at which the system hazard rate reaches its 


maximum, Te can be found by solving 
A, e + A. e - (A, - A ¢ came 
a i I. 2 


for t. While this can be solved numerically, a useful 
analytic sqimmntion if Oia Ts has not been found. However, it 
was discovered that the time at which the system hazard 
fase recaches itS maximum 1S approximately related to the 


time tT at which Po (t) reaches its maximum by 


hronm Che nNumerTeaqiareswlvs el taappcars Umar Ts is 
related to tT by some function of a = ~ ; ene imbalance 
i 
between the two components, i.e. T. = SCO) 8 )t.  Thetexater 


form of g(a) is unknown but it appears to be some function 
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whose values are plotted in Figures 2, 3, and 4. Note 
that g(a) is never greater than 2.17, nor less than 2.0 


fomewall 1.000 << eec 00) yande tor varies ol a > 15; 
ere < 2.01. 
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